Electronic structure and magnetic properties of Li 2 ZrCu04 - a spin 1/2 Heisenberg 

system in vicinity to a quantum critical point 
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Based on density functional calculations, we present a detailed theoretical study of the elec- 
tronic structure and the magnetic properties of the quasi-one dimensional chain cuprate L^ZrCuCU 
(Li2CuZr04). For the relevant ratio of the next-nearest neighbor exchange J2 to the nearest neigh- 
bor exchange J\ we find a = —J2/J1 = 0.22 ± 0.02 which is very close to the critical point at 
1/4. Owing this vicinity to a ferromagnetic-helical critical point, we study in detail the influence 
of structural peculiarities such as the reported Li disorder and the non-planar chain geometry on 
the magnetic interactions combining the results of LDA based tight-binding models with LDA+J7 
derived exchange parameters. Our investigation is complemented by an exact diagonalization study 
of a multi-band Hubbard model for finite clusters predicting a strong temperature dependence of 
the optical conductivity for Li2ZrCu04. 
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I. INTRODUCTION 

Low dimensional magnetism has always attracted great 
interest in solid state physics and chemistry. Especially, 
the phase diagram of low dimensional spin-1/2 lattices 
has been investigated extensively, both theoretically and 
experimentally. For low dimensions and spin-1/2, the in- 
fluence of quantum fluctuations becomes crucial for the 
ground state of the system. The role of quantum fluc- 
tuations is even more pronounced if the system under 
consideration exhibits strongly frustrated interactions. 

Taking only nearest neighbor interactions into ac- 
count, pure geometrical frustration in two dimensions 
(2D) can be realized by special symmetries, i.e trian- 
gular or Kagome lattices. Prominent real material re- 
alizations for such lattices are k-(BEDT-TTF) 2 X (RefQ) 
and ZnCu3(OH)gCl2>2, Another possibility to realize frus- 
trated couplings arc competing nearest neighbor (NN) 
and next nearest neighbor (NNN) interactions. This way, 
the unfrustratcd NN Heisenberg square lattice becomes 
frustrated by adding antiferromagnetic (AFM) interac- 
tions to NNN. The ground state of this model is deter- 
mined by the ratio a=— J2/ J\ between NNN and NN ex- 
change interaction. Including a tiny, but non-vanishing 
coupling between the magnetic planes that is always 
present in real materials, for small a an antiferromagnetic 
ground state is observed like in all undoped HTSC parent 
compounds.— For large a the so called columnar order is 
expected, as found in Li 2 VSi0 4 or Li 2 VGe0 4 ^^ Be- 
tween these commensurably ordered phases, theory pre- 
dicts a spin liquid ground state. The possible realization 
of this spin liquid ground state in PbVC>3 is presently 
under debate 

In one dimension (ID), pure geometrical frustration 
due to NN exchange is impossible, but frustration by 
competing NN and NNN exchanges may occur in close 
analogy to the 2D square J1-J2 model. The phase dia- 
gram of this seemingly simple model is very rich. De- 



pending on the frustration ratio a, a variety of ground 
states was observed in corresponding quasi-lD systems: 

(i) ferromagnetically ordered chains in Li 2 Cu02^ (0 < 
a < 0.25, within an effective single-chain analysis ig- 
noring the antiferromagnetic intcr-chain interaction 12 ), 

(ii) helical order with different pitch angles in L1VCUO4, 
LiCu 2 2 , and NaCu 2 02±2»I£i±£iI&iLI£ (a > 0.25) or spin 
gap behavior like in the famous spin-Peierls compound 
CuGeO^ (—a > 0.241), i.e. both exchange couplings 
are antiferromagnetic. 

Very recently, large interest focused on chain systems 
that are close to the quantum critical point (QCP) at 
a = 0.25. In the vicinity of a QCP, the system is expected 
to answer in a pronounced way to small external per- 
turbations like magnetic fields or pressure. In addition, 
even very small additional exchange paths, especially the 
inter-chain coupling, become important for the magnetic 
properties of the system. A most promising compound 
to study the close vicinity to the QCP is Li 2 ZrCu04, 
where a frustration ratio of a = 0.29 was reported from 
the evaluation of susceptibility and specific heat data^ 
The expected strong field dependence of the specific heat 
from magnetic fields up to 9 Tesla was observed, although 
the NN exchange J\ was estimated to be rather large in 
the order of 300 K. The vicinity of Li 2 ZrCu04 to the 
QCP was supported by preliminary band structure cal- 
culations that obtained a frustration ratio a ~ 0.23. 

In this paper, we report a detailed electronic structure 
study based on density functional calculations. In par- 
ticular, we investigate the influence of the experimentally 
observed Li disorder— for one of the crystallographic Li 
sites on the magnetic properties. We discuss the depen- 
dency of the exchange integrals and the frustration ra- 
tio a on the strong Coulomb repulsion within the Cu 
3d orbitals. In addition, we elucidate the crucial impor- 
tance of the structural distortion of the Cu0 2 chains in 
Li 2 ZrCu0 4 . This distortion is responsible for a reduction 
of both the inter-chain coupling and the NNN exchange 



J 2 , as a consequence placing the system very close to the 
QCP. 



II. METHODS 

For the electronic structure calculation a full- 
potential non-orthogonal local-orbital minimum-basis 
scheme (FPLO5.00-19 and FPLO7.00-28),-22i23 within the 
local (spin) density approximation (L(S)DA) was used. 
In the scalar relativistic calculations the exchange and 
correlation potential of Perdew and Wang was used.— We 
also applied the general gradient approximation (Pcrdew- 
Burkc-Ernzcrho^,) for the exchange and correlation po- 
tential to check whether this influences the LDA re- 
sults. Neither for the band structure nor for the total 
energy differences we found any significant changes. For 
the basis the following states were taken into account: 
Cu (3s3p)/4s4p3d, Li ls/(2s2p) + 3d, O (2s2p3d), Zr 
(As^p) / (bsbpAd) (notation: semi-core states / valence 
states). The Cu 3s3p, Li Is and Zr 4s4p states were 
treated as semi-core states due to the large extension of 
their wave functions. The unoccupied states Li 2p3d, O 
3d and Zr 5p were considered as hybridization states to 
complete the basis. All lower laying states were treated 
as core states. To ensure accuracy of the total energy 
300 k-points in the irreducible part of the Brillouin zone 
were used. 

To model the Li(l) split position of the system (see 
below), we employed the virtual crystal approximation 
(VCA) and coherent potential approximation (CPA). Al- 
ternatively super-cells with different Li(l) orders have 
been studied. To treat the strong on-site Coulomb re- 
pulsion for Cu 3d orbitals explicitly, the L(S)DA + U 
method 2 ^ with U^d m the representative range from 
5.5 eV to 8.0 eV and / = 1 eV for the intra atomic ex- 
change have been used. 

To study some aspects of the magnetic properties of the 
system a tight binding model (TBM) was derived from 
the LDA-calculations to find the transfer integrals for 
a single-band effective extended Hubbard model which 
subsequently has been mapped onto a corresponding 3D 
Heisenberg model. 

In addition, exact diagonalization studies of a five- 
band pd extended Hubbard model and its mapping onto 
the lD-Hcisenberg model using finite open Cu„02n+2 
clusters have been carried out. 



III. STRUCTURE 

The 7-phase of Li2ZrCuC> 4 21 ' 27 crystallizes in an or- 
thorhombic space group and is shown in Fig. QJi. The 
compound contains distorted CuC>4 plaquettes which are 
arranged as edge-shared CuC>2 chains with an Cu-O-Cu 
bond angle of 94.13°. These edge shared CuC>2 chains 
form layers together with the Li(2) atoms in the yz-planc 
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FIG. 1: (Color online) (a) Crystal structure of Li2ZrCu04: 
Chains of edge-shared CUO4 plaquettes (yellow) run along 
the z-direction and are connected by ZrOe-octraheda (light 
blue). The fully occupied Li(2) site (gray) lies in the chain 
plane. The Li(l) split position (pink) alternates with the 
ZrOe octahedra. (b) Section of the Cu-0 chain: The O atoms 
of a plaquette are deflected from the mirror plane (shown in 
yellow) . 

interconnected by ZrOg octrahedra. The alternating ar- 
rangement of the connecting ZrOe octrahedra lead to the 
deviation of the Cu02 chains from the ideal planar ge- 
ometry. In Fig. [TJa the deflection of the O atoms from 
the chain plane is illustrated in detail. 

As the CuC>4 plaquettes exhibit a localized, effective 
spin 1/2, their linking and surrounding is of main im- 
portance for the formation of the magnetic ground state. 
A first structural analysis suggests sizable ferromagnetic 
NN exchanges caused by the close to 90° Cu-O-Cu bond 
angle according to the Goodenough-Kanamori- Anderson 
rule, provided the direct Cu-Cu transfer integral can be 
ignored.— 

A further structural characteristic of the 7-phase is 
a split position for Li(l) which is placed between the 
chain layers. The distance of the Li(l) atoms from the 
high-symmetry position^ is 0.37A. While the influence 
of the split position of Li(l) (corresponding to 50% disor- 
der within a classical picture, ignoring possible tunneling 
processes of Li(l) between the two sites of the split posi- 
tion) to the magnetic properties is investigated carefully, 
we neglect the small disorder of 3-5 % at the Li(2) and 
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Cu site in our theoretical calculations as a good approx- SYM 
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It is well known that in the vicinity of a QCP small 
changes in any parameter can be of crucial relevance for \ / \ / / \ 

the ground state realized by the system. Therefore, the 
above mentioned structural peculiarities of I^ZrCuOg 

raise the following two questions especially related to ^^..^--^=^-zrJ^--^-T^~^^^p^^-T^~~^% 
their influence on the electronic and magnetic proper- \ / \ / \ / \ / \ 

ties. (i)How does the modeling of the Li(l) split position 
influences the calculated ground state? As there exists 
no standard procedure in band structure codes to treat 
split positions we suggest various classical (i.e. within the 
adiabatic approximation) approaches, (ii) To which ex- 
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the placement of the compound in the J\ — J<z phase dia- 
gram? Therefore a fictitious structure with ideal planar / | / \ ./ \ / / \ / 

chain geometry was also considered^ Qy Q 



IV. RESULTS AND DISCUSSION 

A. Modeling the Li(l) split position 

The vicinity of Li2ZrCu04 to a QCP requires a careful 
check of the structural description of the compound as a 
basis for our theoretical calculation. Therefore the treat- 
ment of the Li(l) split position might be crucial, although 
the different models may only exhibit subtle differences. 
In the following we suggest different models to handle 
the split position and discuss the results with respect to 
the influence of these model assumptions on the relevant 
states and related dispersions. 

(i) To start from the most simplest approximation 
we performed electronic structure calculation for Li(l) 
placed at the high-symmetry position (see Fig. [5] (SYM)). 

(ii) As an alternative structural model we carried out 
VCA calculations (see Fig. [3] (VCA)). This approach oc- 
cupies both sites of the split position simultaneously by 
a half of a Li atom. Each of this modified Li atoms is 
constructed by carrying the half number of valence elec- 
trons and a modified core to conserve charge. For the 
VCA calculations the distance between the two sites of 
the split position has to be enlarged to suppress the over- 
lap of the Li valence states and ensure the convergence of 
the calculation^ Also asymmetric charge distributions 
between the two sites of a split position (as e.g. 0.8 and 
0.2 Li) lead to nearly identical dispersion of the relevant 
states as long as the total charge of one Li atom per split 
position is conserved. However, deviations of the total 
charge of a split position cause sizable changes in the 
electronic structure. 

This is also the reason why the description of the Li(l) 
split position by CPA fails, caused by the loss of charge 
balance at the split position in the later approximation. 
In contrast to the VCA, where the total charge at each 
split position sums up to one Li atom, the random occu- 
pation of the Li(I) split sites in the CPA also generates 
cases with empty and double occupied Li(I) sites with 
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FIG. 2: (Color online) Different structural models to treat 
the Li(l) split position. From top to bottom: SYM - high- 
symmetry position for Li(l), VCA - both sites occupied by 
partial Li(l) atoms, SI to S3 - super-cells with different Li(l) 
pattern. In SI and S2 the Li(l) atoms along the chains are 
shifted parallel while in S3 the Li(l) atoms alter along the 
chain direction. SI and S2 differ in the Li(l) inter-chain ar- 
rangement. The slight tilting of the ZrOe-octrahedra (see 
Fig. [TJ is not depicted for the sake of simplicity (full lines 
foreground, dashed lines background). 
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FIG. 3: (Color online) Calculated LDA density of states and 
band structure for the physical relevant states using different 
structural models for the Li split position. The comparison 
shows just minor differences. 



dramatic influence on the band structure and the total 
energy. 

(iii) As both latter approaches do not touch the local 
symmetry, within a further model super-cells with dif- 
ferent Li(l) patterns have been calculated, breaking the 
local symmetry at the split position but suffer from long 
range order (see Fig. [2] S1-S3). In the two upper pan- 
els the super-cells were constructed by a reduction of the 
space group symmetry. The Li(l) atoms along the chains 
arc all shifted in the same direction while the Li(l) atoms 
of neighboring chains are parallel (SI) or anti-parallel 
(S2) shifted. The third super-cell (Fig. [2] S3) shows an 
alternating displacement of the Li(l) atoms along the 
chain obtained by a doubling of the cell in y-direction. 

The resulting bands and density of states for the differ- 
ent structural models show a very similar total behavior 
(see Fig. [3]) . Small deviations in VCA calculations can 
be observed in the bonding region of the valence band 
around -5.5 and -2.0 eV (not shown) in comparison to 
the high symmetry model. These differences can be un- 
derstood by the influence of the two modified Li cores to 
the crystal potential and the corresponding shift of Li-0 
states. 

Comparing the density of states and LDA band struc- 
ture for the anti-bonding dpa* states, which are relevant 
for the magnetic exchanges and therefore the magnetism 
of the system, we found only tiny differences, (i) The 
VCA model results in an slightly larger bandwidth com- 
pared to the other models, (ii) The high symmetry and 
super-cell models nearly fit perfect. The additional bands 



of the S3 super-cell with an alternating Li pattern can be 
understood due to the doubling of the primitive cell hav- 
ing its origin in the symmetry breaking of the Li pattern. 



B. Electronic Structure and Magnetic Properties 

As the different structural models for the Li(l) split 
position show no significant difference in behavior for 
the relevant low energy states (see Fig. [3]), the electronic 
structure and the microscopic magnetic model are dis- 
cussed for the SYM model as a good representative for 
all different structural models. To estimate more quanti- 
tatively the differences caused by the choice of the struc- 
tural model the results for the SYM and the VCA models 
are compared. 
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FIG. 4: (Color online) Calculated atom and orbital resolved 
density of states (LDA) for Li2ZrCu04 applying the SYM 
model for the Li(l) split position (see text). The Li contri- 
bution to the valence band is negligible. Cu and O dominate 
the anti-bonding dpa* band. 

In Fig.Q]the atom and orbital resolved density of states 
for Li2ZrCu04 are shown. The valence band has a width 
of about 7eV, comparable with other chain cuprates and 
is dominated by Cu and O states, especially at the lower 
edge of the valence band around -6 eV (bonding states) 
and at the energy range close to zero (anti-bonding 
states). States of the Zr-0 octrahedra appear in the mid- 
dle part of the valence band. Their contribution to the 
states close to the Fermi energy is negligible. 

The metallic character of the LDA results is in contra- 
diction to the experimentally observed insulating ground 
state. This effect is a well known shortcoming of the LDA 
as it underestimates the strong correlations in the Cu 
3d states. The related strong Coulomb repulsion on Cu 
sites can be considered explicitly by a model approach 
(e.g. a single-band Hubbard model using tight-binding 
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FIG. 5: (Color online) Valence band together with the band 
characters of the Cu 3d and O 2p anti-bonding states (blue) 
indicated by the line width. The band at Fermi energy is 
dominated by Cu and O anti-bonding dpa* states. 



FIG. 6: (Color online) Left: Comparison of the calculated 
band structure (black dots) and the least-square fit (red line) 
of the TB model. Right: Considered hoppings for the TB 
model. The leading interactions are labeled and highlighted 
by red colors. 



parameters derived from the LDA) or by LDA+J7 cal- 
culations. Both methods were applied in this study and 
lead to an insulating ground state in agreement with the 
experimental data. Nevertheless, the LDA calculation 
delivers important insights into the electronic structure 
and provides the relevant bands for the development of 
a microscopically based picture for the magnetic ground 
state of Li 2 ZrCu04. 

The well separated states around the Fermi energy are 
of main interest as they determine the magnetism of the 
system. These states arc dominated by the in-plane or- 
bitals of the CUO4 plaquettes corresponding to the anti- 
bonding dpa* states (see two lower panels in Fig.|4]). The 
distortion of the plaquettes from an ideal planar chain ge- 
ometry does not affect significantly the band characters 
of the anti-bonding dpa* states. The hybridization with 
out-of-plane states is very small and does not differ from 
compounds with ideal planar Cu02 chains. 

In Fig. [5] the band structure of the valence band to- 
gether with the Cu 3d x 2 _ y i and O 2p x ,y characters are 
shown. The anti-bonding dpa* bands (see Fig. [5] left) 
show their main dispersion along the chain direction (Z- 
r) and in the chain plane (r-Y, S-X) indicating the lead- 
ing magnetic exchanges in the layers. Perpendicular to 
these directions the dispersion is small, suggesting small 
exchanges between the chain layers. 

To unravel the hierarchy of exchanges in the system 
as well as to estimate a cut-off for negligible long range 
exchanges in further super-cell calculations, an effec- 
tive one-band model was fitted to the anti-bonding dpa* 
bands. Effective hoppings for Cu-Cu distances up to 12 A 
are considered (see Fig. |5] right) and calculated using a 
least-square fit procedure. The resulting TB fit and the 
LDA bands are in perfect agreement (see Fig. |S| left). 

All obtained transfer integrals larger than 5 meV are 
summarized in Table [IJ We found the main hoppings be- 
tween NN and NNN along the chains and three consider- 
able inter-chain hoppings in the chain plane. Comparing 
the SYM results with results from the VCA model, we 
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TABLE I: Resulting hopping integrals from TB fits for differ- 
ent structural models. 



see that the differences in the hopping integrals are very 
small and in the same order as the accuracy of the TB 
fit procedure. This indicates that the chosen structural 
model for the Li(l) split position is irrelevant for the cal- 
culated coupling constants of Li2CuZr04. Taking into ac- 
count the strong Coulomb repulsion in a model approach 
we construct from our TB model a single-band Hubbard 
model (with a typical valued U e ff = 4.2 eV for edge 
shared CuC>2 chains) which can be mapped subsequently 
onto a Heisenberg model. We find J^ FM = 1.7±0.1meV 
for the NN exchange and J 2 4FM = 3.0 ± 0.1 meV for the 
NNN exchanged 

The exchange parameters resulted from the effective 
TB model disregard the FM contributions to the ex- 
changes which are, however, expected between NN pla- 
quettes in the compound due to the close to 90° bond 
angle Cu-O-Cu along the chains. Naturally, these FM 
interactions are included in LDA+[/ calculations. To 
evaluate J\ and J2 on a quantitative level, we performed 
a series of super-cell calculations. The comparison of 
the total energies of these super-cells with different spin 
arrangements^ and a subsequent mapping onto a Heisen- 
berg model leads to a FM NN exchange J™ T = -11.2 
meV and an AFM NNN exchange jJ OT = 2.2 meV. The 
ground states of the compound is determined by the ra- 
tio a = -jJ OT /jf OT = 0.22 ± 0.03 that would re- 
sult in a FM ordered ground state in close vicinity of 
a QCP (a = 0.25). The variation of U 3d from 5.5 eV 
to 8.0 eV in the LDA+U calculation does not change 
this behavior qualitatively, although for the smaller U 
values a slightly increased a is obtained^ Besides a 
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TABLE II: Total and partial exchange integrals for different 
structural models from TB model and LDA+C/ calculation 
(U 3d = 7.5eV). 
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FIG. 7: (Color online) Comparison of band structure and 
density of states at Fermi energy for both chain geometries. 
The small displacement of O atoms from the mirror plane 
leads to a decrease of the band dispersion. 



possible small error in the calculated exchange parame- 
ters, the experimentally observed spiral state from NMR 
data— might result due to the sizable inter chain cou- 
plings J$ = A{tl) 2 /U cS ~ 0.6 meV and Jf - 0.6 meV 
that both stabilize long range order. 



shown. For a planar chain geometry we obtain an about 
50% larger bandwidth, but a very similar shape of the 
bands compared with the real structure. On first glance, 
this suggests a scaling of the leading hopping parameters, 
only. Our TB fit yields nearly a doubling of the hopping 
integrals as listed in Tab. U thus a dramatic change of 
the corresponding exchange terms can be expected. To 
take into account the large FM contributions to J™ T , 
we carried out LDA+C/ calculations for the same super- 
cells used before— 

Surprisingly, we find a nearly unchanged NN exchange 
J TOT , but only due to the compensation of the in- 
creased individual contributions J[ M and J^ FM (see 
Tab. HH. Since the FM contribution J 2 FA/ to jJ OT re- 
mained small - as for the real structure - it basically 
scales like J^ FM ~ ^/Ucs- Thus, we obtain a large 
change in jJ OT and consequently in the frustration ra- 
tio a ~ 0.6 for the planar chains. In consequence, the 
ratio a = —J2/J1 depends strongly on structural details 
of the local Cu-0 environment. This is especially impor- 
tant regarding the vicinity of the system to the QCP. On 
the other hand this sensitivity of a to the chain buckling 
may provide the opportunity of manipulating the ground 
state of the system selectively by substitution or external 
pressure. 

Unfortunately, a planar chain arrangement leads also 
to a sizable increase of the inter chain coupling. This 
can be understood as the decrease of buckling increases 
the inter-chain overlap of the O orbitals belonging to the 
pda* states. This way, a tendency towards long range 
order is stabilized, although the quantitative influence 
of inter-chain couplings to the ground state is sparsely 
considered in the literature. 



D. Stability of the Li(l) split position 



C. Influence of chain-buckling 

Whereas the different treatment of the Li(l) split po- 
sition does not influence the band structure and the cor- 
responding magnetic model essentially, the role of the 
distortion of the CUO4 plaquettes (see Fig. [TJd) remains 
to be addressed: How much does the non-planar geome- 
try of the chains influence the balance of the NN, NNN 
and inter-chain couplings and the resulting placement 
of the compound in the J\-J% phase diagram, especially 
its vicinity to the QCP? For a fictitious planar chain^ 
changes are expected since a changed Cu-O-Cu bond an- 
gle (changed by about 0.8°) will influence J TOT and the 
different orientation of the anti-bonding pda* "molecular 
orbital" will alter J™ T . Since these changes are hardly 
predictable quantitatively from structural considerations 
only, we apply our mapping approach also to a fictitious 
compound with ideal planar Cu02 chains. 

In Fig. [7]thc resulting anti-bonding dper* bands and the 
related density of states for this fictitious structure are 



Although the different treatments of the Li(l) split po- 
sition yield no significant influence on the magnetic ex- 
change parameters the split position may influence other 
properties of the compound. Related to the split posi- 
tion, the local symmetry or disorder could modulate the 
pitch angle of the helical state or be crucial for the for- 
mation of a possible multiferroic phase. Furthermore, 
7 Li NMR and complex dielectric measurements reveal a 
glass-like ordering of the Li(l) ions below T g = 100 K. 37 
Thus, using our super-cells simulating different types of 
Li(l) order (see Fig. [2]), we tackled the questions whether 
a tendency to a static Li(l) order can be supported by 
calculations. 

For the purpose of computational feasibility, we re- 
stricted ourselves to a fourfold cell (32 atoms).— The 
comparison of the total energy for the different ordered 
Li(l) patterns (SYM, SI, S2, S3) favors energetically the 
high-symmetry model (SYM) with a single Li(l) posi- 
tion. The second lowest in energy is the structure with 
alternating Li(l) displacement along the chain direction 
(S3). For this Li(l) arrangement, the dependence of the 
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displacement of Li(1) from high-symmetry position (dx/A) 

FIG. 8: (Color online) Energy difference depending on the 
displacement of the Li(l) atom from the high-symmetry po- 
sition. The experimental observed displacement of 0.4 A is 
indicated by an arrow. The inset shows a harmonic (second 
order) fit to the calculated energy differences for small dis- 
placements from the equilibrium position (±0.lA). The same 
fit is shown in the main panel, indicating deviations from the 
harmonic behavior at larger Li(l) elongations. 



total energy on the displacement is depicted in Fig. [5] 
The resulting curve does not indicate a double well po- 
tential as expected for a Li(l) split position, but a mini- 
mum around the high symmetry position and a harmonic 
behavior (see inset Fig. [5]) up to almost 0.3 A Li(l) dis- 
placement. According to our calculations, for the ex- 
perimentally suggested Li(l) split position (marked by 
the arrow in Fig. [5]) only slight anharmonic effects could 
be expected. We also investigated whether the strong 
Coulomb repulsion on the Cu site could be the reason 
for the discrepancy between the single Li(l) position sug- 
gested by the calculations and the experimental split po- 
sition using LDA+C/. With U^d in the range from 5 - 
7.5 eV we do not obtain any significant changes of the 
position of the minimum and the slope of the curve, sug- 
gesting that the correlations at the Cu site play only a 
minor role for a possible Li(l) order. 

So far our calculations with a fixed lattice (apart from 
Li(l)) can not support the split position, but favor the 
high-symmetry solution (SYM). The discrepancy with 
the experiment that observed a glass-like Li(l) ordering 
below 100 K, may have the following origins: (i) dynamic 
effects that violate the adiabatic approximation used in 
the calculations or (ii) a correlated Li(l)-0(1) position 
without long range order that would be difficult to de- 
tect in X-ray diffraction experiments. An investigation of 
these effects is beyond the scope of the present study but 
should be carried out in future since it might be crucial 
for the final picture of the compound. 



E. Cluster Calculations 

In order to provide an independent microscopic in- 
terpretation of the exchange integrals reported above, 
we performed exact diagonalization studies for the com- 
monly used five-band extended pd Hubbard model for 
finite open Cu„02n+2 clusters. The most reliable way 
to choose a reasonable parameter set by fitting various 
spectroscopic experiments cannot be used here due to the 
lack of single crystals and available experimental data. 
Since no parameter set for the five-band pd model for a 
buckled chain geometry is available, we considered only 
the fictitious planar case (FIC) without any O-buckling. 
This way, a direct comparison to the band structure de- 
rived results (FIC) should be more reliable. Thus, we 
started from an available parameter set known for the 
well-studied closely related sister compound Li2Cu02 
with planar Cu02 chains i 39 i 40 Then, only two Hamil- 
tonian parameters - the mean Cu-0 on-site difference 
A pd = 3.5 cV (3.75 eV) and K pd = 54 meV (81 meV) 
- have to be slightly lowered to obtain consistency with 
the band structure results for the exchange integrals J\ 
and J2 reported above. For comparison the Li2Cu02 pa- 
rameters are given in brackets^. The slightly smaller 
Apd is in line with the slight increase of the Cu-0 dis- 
tance by about 2% along the Cu02 chains for Li2ZrCu04 
compared to Li2Cu02- 

With these parameters the title compound exhibits a 
charge gap of about 2.35 eV as expected for an undopcd 
cuprate. Therefore all eigenstates below that energy are 
spin states which can be described also reasonably well 
by a Heisenberg Hamiltonian. In case of a CU5O10 cluster 
we deal with 10 spin states. Mapping for instance these 
10 spin states onto a ID Ji- J 2- J3-Heisenberg model^ 
we arrive then at J\ = —10.6 meV, J2 = 5.9 meV and 
J3 = —0.3 meV close to the values (FIC) given in Tab. II. 
Notice that our mapping on larger clusters with 5 Cu- 
sites is considered to be more accurate than the determi- 
nation of J\ from a Cu204-dimer and J2 from a CU2O8 
pscudo-trimer, where the central Cu has been removed^ 

In the context of the lower value of J\ compared with 
Ref. [2^ we would like to mention that from a micro- 
scopic point of view J\ depends strongly on the adopted 
value of the direct ferromagnetic Cu-0 exchange param- 
eter K p d- Within perturbation theory it occurs already 
in the second order whereas the antiferromagnetic su- 
perexchange due to deviations from the 90° and/or due 
to different crystal fields seen by the O orbitals along 
and perpendicular to the chain direction as well as the 
ferromagnetic Hund's-rule coupling and the O-sites oc- 
cur in the fourth order, only4^ Analytically and numer- 
ically one obtains a linear dependence. Unfortunately, 
at present there is no consensus about a reliable value 
of Kpd based on microscopic considerations. Therefore 
it is frequently treated as fitting parameter to reproduce 
empirical or otherwise determined values of J\ just as 
done here. For other Cu02 edge-shared chain cuprates 
.Kpd-values in between 50 meV and 110 meV have been 
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FIG. 9: (Color online) Optical conductivity of Li2ZrCu04 
within the five-band extended Hubbard model for a CU5O12 
cluster. 



employed by various authorsi 20 i 39 i °i 3 i 5 In our previous 
work on Li 2 ZrCuO,f2S 70 meV has been adopted, thus 
with the given above 54 meV we are near to the lower 
bound of so far used -Kp^-values. The former value was 
based on ID-fits of the magnetic susceptibility x(T) and 
of the magnetic specific heat c v (T) ss c p (T), which might 
be significantly affected by the ignored inter-chain in- 
teraction. Due to the vicinity to the critical point the 
calculated ID saturation field at T = OK is very small 
of the order of 4 to 5 Tesla. Hence, the experimentally 
measured saturation field of about 30 Tesla (at low tem- 
perature) is clearly dominated by the inter-chain inter- 
action in accord with our calculated sizable inter-chain 
couplings (see Sec. lIVBl) . 

In order to stimulate optical measurements and for a 
further experimental check of the applied Hamiltonian 
parameters we calculated (see Fig. the in-chain opti- 
cal conductivity for 300 K and 7 K using the technique 
presented in Ref. H^. We draw attention to the strong 
temperature effect near 3 eV between low (7 K) and room 
temperature (300 K). This phenomenon is related to the 
different optical response in different initial spin states. 
For relatively small exchange integrals J\ and J2 the ther- 
mal population of these spin states changes markedly 
in the mentioned temperature range between zero and 
room temperature. Since the NN exchange J\ (FIC) 
is mostly unchanged compared to the buckled structure 
(SYM, VCA), a strong temperature dependence of the 
optical conductivity should hold for the real compound. 

Further investigations to settle the value of 3\ are 
highly desirable not only to get access to the Hamiltonian 
parameters,— but also to check for renormalization and 
creation of longer-ranged couplings due to strong enough 
spin-phonon interaction^ The strong coupling of the ex- 
change parameter J2 to the O buckling, indicated by our 



band structure calculations (sec Tab. HI]) , opens room for 
such a scenario. 



V. CONCLUSIONS AND SUMMARY 

To summarize, applying LDA and LDA+J7 band struc- 
ture calculations as well as exact diagonalization studies 
of an extended pd Hubbard model for large Cu n 02 n +2 
clusters (n=5), that were finally all mapped onto a 
Heisenberg model, we obtain a consistent, microscopi- 
cally based picture: (i) Li2ZrCu04 can be understood as 
a quasi-lD chain compound with FM NN J\ and AFM 
NNN J2 exchange in close vicinity to a QCP. (ii) We find 
sizable inter-chain couplings that should be relevant for 
the magnetic ground state of the system, especially due 
to the vicinity of the system to a critical point where 
the impact of the leading exchange parameters on the 
ground state gets small, (iii) Calculations for a fictitious 
structure with planar chains indicate that the balance 
between J\ and J2 depends heavily on structural details, 
especially on the magnitude of the Cu02 chain buckling, 
(iv) This also applies for the inter-chain coupling which 
is strongly increased for the fictitious planar geometry. 

Our findings suggest that the buckling of the CuC>2 
chain is crucial for the vicinity to the QCP, therefore it 
should be possible to tune the system towards that point 
by chemical substitution at the Zr site or by external 
pressure. On the other hand, it would be desirable to 
decrease the inter-chain coupling for which the orienta- 
tion of the CuO pda orbitals is more relevant than inter- 
atomic distances according to our results. Thus, a search 
for good quasi ID systems should turn towards crystal 
structures where the chains are arranged in a strongly 
non-planar pattern. 

Our simulations indicate that the Li(l) split position 
has only a negligible direct influence on the spin system, 
assuming that all other atoms of the structure are fixed 
to their experimental position. In particular, our calcu- 
lations do not support the experimental observation of a 
glassy-like order for Li(l) on this split positional since it 
yields a single Li(l) position as most energetically favor- 
able. However, an indirect influence of the Li(l) position 
on the spin system is possible: Taking into account a 
possible relaxation of the neighboring oxygen atom in the 
Cu02 chain could result in a sizable modification of the 
exchange parameters as shown by the strong influence of 
this O position on the leading exchange terms. Another 
possible origin of the discrepancy concerning the Li(l) 
position could arise from dynamic and non-adiabatic ef- 
fects that are beyond the scope of the present study. 
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